Melting temperature of screened Wigner crystal on helium films by molecular 

dynamics 
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Using molecular dynamics (MD) simulation, we have calculated the melting temperature of two- 
dimensional electron systems on 24uA-500A helium films supported by substrates of dielectric con- 
stants e s = 2.2 — 11.9 at areal densities n varying from 3 x 10 9 cm -2 to 1.3 x 10 10 cm -2 . Our results 
are in good agreement with the available theoretical and experimental results. 
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At suficiently low densities and temperatures, an elec- 
tron gas is expected to undergo a phase transition to a 
lattice (because of the domination of the Coulomb in- 
teraction energy over the kinetic energy) which has re- 
ceived the name Wigner crystal^. The two-dimensional 
(2d) Wigner crystal is well established and experimen- 
tally it was first observed on liquid helium surface |2j and 
more recently in semiconductors structures like MOS- 
FET's and heterojunctionsQ . These systems can be used 
for testing several theoretical predictions in many-body 
theory, such as phase-transitions of the electron system, 
metal-insulator transition and now electrons on helium 
surface are being proposed as a set of strongly interact- 
ing quantum bits for quantum computers Q. Electrons 
on the surface of bulk helium form a crystal at a tem- 
perature T m = 2e 2 (7m) 1 / 2 / '{ejie + l)r m , which is much 
higher than the Fermi temperature, Tp — nnh 2 jra in a 
density range of 10 5 — 10 9 cm -2 (where n is the elec- 
tron areal density, ejy e is the dieletric constant of helium 
and r m is the plasma parameter in the melting temper- 
ature defined as the ratio of potential to kinetic energy) . 
Therefore, such electrons in this regime obey the clas- 
sical Boltzmann statistics. Experimentally the liquid to 
solid transition in the bulk takes place for a value of the 
coupling constant T m = 137 ± 150 and computer sim- 
ulations of Kalia et al. [f| showed an agreement with the 
experimental measurements indicating a first-order melt- 
ing in T m = 118 - 130. 

Superficial electrons on liquid helium films form also 
a very interesting system to study the many-body prop- 
erties of 2d screened systems. In this case the screen- 
ing is provided by the image charges in the substrate 
beneath the film. The screening effect can drastically 
change the electron-electron interacting potential going 
from 1/r to 1/r 3 through varying external parameters 
such as film thickness and dielectric constant of the sub- 
strate. Peeters using a phenomenological approach 
got a reduction in the phase diagram of this electron sys- 
tem comparing with the bulk case. SaitohQ obtained the 
melting transition in this system using an analytical ap- 
proximation to the angular frequency of the transverse 
Wigner phonon combined with the Kosterlitz-Thoules 
melting criterion. His result is in agreement with the 



experiment by Jiang et a/.[j|. Candido et al. 8] studied 
the thermodynamical, structural and dynamical prop- 
erties of this two-dimensional electron system by com- 
puter simulation. Experimmentally, the melting temper- 
ature of the Wigner crystal on thin helium films adsorbed 
on dielectric substrates was measured by Jian g et a/.[|| 
through electron mobility and by Mistura et al. 10] using 
microwave cavity technique. 

In this paper, we present an accurate MD calculation 
for the melting temperature for an electron system over 
a helium film adsorbed on a dielectric substrate. In Fig. 
1 we show schematically the geometrical arrangement of 
the system considered. The obtained results are directly 
compared with the available experimental data, of Mis- 
tura et aL[l(j and Jiang et al. 9], and the theoretical re- 
sults of Peeters^ and Saitoh|J 

We consider a two-dimensional system of electrons on 
a helium film of thickness d adsorbed on a substrate 
of dielectric constant e, s , interacting through a screen- 
ing Coulomb potential|ll|. The electron system is im- 
mersed in a rigid, uniform, positively charged background 
to make a neutral charged system. The Hamiltonian for 
such a system is given by 
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where <5 = (e s — l)/(e s + 1) with the dielectric constant 
of helium approximated by 1 (eiie = 1-057) and C/f, is 
the interaction of electrons with the uniform positively 
charged background. 

In this work most of the molecular dynamics calcula- 
tions were performed on a system of 100 electrons with 
a few runs of 484 and 784 electrons to study size ef- 
fect. The finite size effect is investigated by changing 
the system size and the thermodynamical behavior in 
an infinite system is derived from their extrapolation. 
The initial position of the electrons is a triangular lattice 
which is accomodated in a rectangular box with periodic 
boundary condictions to eliminate the surface effects. Be- 
cause of the long range nature of the electron-electron 
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FIG. 1: Schematic view of the electron system. 



and electron-background interacting potential we are em- 
ploying Ewald summation which splits the potential into 
a long-range and a short-range part. The long-range part 
is handled in the k-space, while the short range part in 
the real space. We have used the fifth-order predictor- 
corrector algoritm to integrate the Newton's equation of 
motion with the MD time step varying from 10" 12 to 
10" 15 sec, since it has some scale dependence on the 
electron densities. The optimum time step leads to a 
conservation of the total energy of 1 part in 10 4 after 
several thousand time step runs. The time average of 
the physical quantities were obtained over 120 000 time 
steps after the system has reached equilibrium. 

In Fig. 2 we present the results for the total energy 
per electron versus temperature to illustrate the general 
feature of the melting transition in this system. The solid 
squares in the figure represent the results for the electron 
liquid which has been monotonically cooled from a higher 
temperature. The open circles are the results for the elec- 
tron solid which has been monotonically heated from a 
lattice at very low temperatures. It means that our sim- 
ulations were performed in cascade, i.e., an equilibrated 
configuration obtained for a given higher (lower) temper- 
ature was used as an input to reach another configuration 
at lower (higher) temperature. As one can see the elec- 
tron system shows hysteresis and latent heat on melting, 
which characterize a first order transition as other 2d 
classical systems. The melting temperature range is 1.83 
K < T < 2.05 K defined from the vertical dashed lines in 
Fig. 2 representing the hysteresis region. Thus, we would 
define the melting temperature, T m , as exactly the mean 
point in the temperature width of the hysteresis AT, i.e, 
T m = (1.94 ± 0.11) K with the error bar giving by half 
the temperature width of the hysteresis. The value of the 
latent heat per particle and the change in the entropy on 
melting are found to be 0.40 K and 0.21 ke, respectively. 
We also find our MD results for the melting temperature 
is in agreement with those of Kalia el aZ.Q for the bulk 
limit. 

Fig. 3 shows size dependences of the transition tem- 
perature T m for different electron densities. The error 
bars on T m indicate the hysteresis width. When elec- 
tron number becomes larger, T m decreases because the 
periodic boundary conditions favors the solid phase. The 
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FIG. 2: Total energy per electron as a function of temperature 
for a system of N = 100 electrons on a helium film supported 
by a glass substrate, e s — 7.3, film thickness d = 240 A and 
density n = 1.3xl0 10 cm -2 . 
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FIG. 3: Size dependence of the melting temperature for elec- 
trons on helium films above a substrate with dielectric con- 
stant t a = 7.3 at four different densities: n = 1.3 x 10 10 cm -2 
and d = 240 A(squares); n = 1.0 X 10 10 cm" 2 and d = 285 
A(triangles-up); n = 0.9 x 10 10 cm" 2 and d = 260 A(triangles- 
down); and n = 0.75 x 10 10 cm" 2 and d = 305 A(diamonds). 
The lines are linear fit. 



transition temperature, however, seems to follow a linear 
decrease as a function of 1/N. Therefore the melting 
point in the thermodynamic limit can be obtained defi- 
nitely by extrapolating the finite size data. 

The extrapolated melting temperature is exhibited in 
Fig. 4 as a function of the electron density (top), film 
thickness (middle) and dielectric constant of the sub- 
strate (bottom). We roughly estimated the error bar on 
the experimental values for the melting temperature to 
indicate the uncertainty of about 15% — 20% on the ex- 
perimental measurement of the electron density. As is 
shown in Fig. 4, our results are in good agreement with 
those obtained experimentally in Ref . [Io| . The top panel 
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FIG. 4: The melting temperature as a function of the electron 
density (top), film thickness (middle) and dielectric constant 
of the substrate (bottom). The experimental results from 
Ref.[To|| are given by solid squares and our MD simulation 
results are indicated by solid triangles. 



trie constant of the substrate leads to a more polarizable 
system with stronger screening. As a consequence, the 
melting temperature goes down. 

In table I, we show a comparison of our MD simulation 
results of the melting temperature with the available the- 
oretical and experimental results for electron systems on 
a thin helium film surface. For densities below 1.0 x 10 10 
cm~ 2 one can see that our MD calculations are in agree- 
ment with both Jiang et al. and Mistura et al. [Tol ] 
experimental measurements. They are also in agreement 
with Saitoh's theoretical results|j| in the range of the 
densities studied. However, we got some discrepancies 
with Peeters's results^ that can be justified, as pointed 
out by Saitoh0, as being a possible double counting on 
T m . For densities larger than 1.0 x 10 10 cm~ 2 , our MD 
simulation melting temperatures are higher than the ex- 
perimental ones, though the differences are almost within 
the uncertainty of the experimental results. A possible 
explanation for this discrepancy is that the quantum ef- 
fect can be important at such densities. Besides, we note 
that the change in entropy on melting decreases as the 
density or the dielectric constant of the substrate (the 
film thichness) increases (decreases). This might imply 
that transition becomes continuous at high densities. 

In conclusion, we have shown that the MD is able to 
reproduce the experimental measurements of the melt- 
ing temperature in the two-dimensional electrons on thin 
liquid He films. Our results are i n g ood agreement with 
those obtained by Mistura et a/. |10| and Jiang et al.\$. 
These results in the classical regime (i.e. n < 1.0 x 10 10 
cm~ 2 ) should be useful to the experimental and theo- 
retical investigation of the melting transition in this sys- 
tem. For larger densities ( n > 1.0 x 10 10 cm -2 and 
e s > 2.2), the results might be beyond the applicability 
of the present method. 



shows that the melting temperature increases with in- 
creasing the electron density (for fixed e s and d) due to 
the the fact that the screening becomes weaker and, con- 
sequently, the electron-electron interaction is enhanced. 
The middle panel also shows a shift in the melting transi- 
tion to higher temperature with increasing the film thick- 
ness (for fixed n and e s ), which is a consequence of de- 
creasing in the screening resulting in a stronger electron- 
electron interaction. In the bottom panel, the melting 
temperature decreases as the dielectric constant of the 
substrate increases at fixed d and n. A larger dielec- 
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